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We present a measurement of the elastic differential cross section da(pp — l pp)/dt as a function 
of the four-momentum-transfer squared t. The data sample corresponds to an integrated luminos- 
ity of « 31 nb" 1 collected with the DO detector using dedicated Tevatron pp Collider operating 
conditions at y/s = 1.96 TeV and covers the range 0.26 < \t\ < 1.2 GeV 2 . For \t\ < 0.6 GeV 2 , 
da/dt is described by an exponential function of the form Ae~ b ^ with a slope parameter 
b = 16.86 ± 0.10 (stat) ± 0.20 (syst) GeV" 2 . A change in slope is observed at \t\ ^ 0.6 GeV 2 , fol- 
lowed by a more gradual \t\ dependence with increasing values of \t\. 

PACS numbers: 13.85.Dz 



I. INTRODUCTION 

The differential cross section, da/dt, for pp — > pp, 
where t is the four- momentum-transfer squared [1], 
contains information about proton structure and non- 
perturbative aspects of pp interactions. In the \t\ range 
studied here, the nuclear scattering amplitude is expected 
to dominate [2], and with increasing \t\, da/dt is expected 
to fall exponentially, followed by a local diffractive mini- 
mum, after which da/dt continues to decrease [3]. 

Studies of da/dt at different center-of-mass energies, 
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y/s, have demonstrated an effect known as shrinkage, 
namely that the slope of the exponential fall-off of the dif- 
ferential cross section as a function of \t\ becomes larger 
with increasing y/s, and the \t\ value at which the local 
diffraction minimum occurs is reduced [4]. It has also 
been observed that the shape of the local diffractive min- 
imum is different between pp and pp elastic scattering [5] . 
The elastic differential cross section plays an important 
role in constraining soft diffractive models [6-8] which 
cannot be directly calculated by pcrturbativc QCD [9]. 

In this Article, we present a measurement of the pp 
elastic differential cross section at y/s = 1.96 TeV in 
the range 0.26 < \t\ < 1.2 GeV 2 , measured using the for- 
ward proton detector (FPD) spectrometer system of the 
DO experiment [10]. Since p and p elastic scattering 
angles are typically very small (on the order of mil- 
liradians), they are not covered by the main DO de- 
tector. The elastically scattered protons and antipro- 
tons are tagged with detectors inserted in the beam 
pipe on either side of the interaction point (IP). Our 
measurement extends the \t\ range previously studied 
by the CDF (0.025 < \t\ < 0.29 GeV 2 ) [11] and by the 
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E710 (0.034 < \t\ < 0.65 GcV 2 ) [12] Collaborations at 
the Tevatron and constitutes the first confirmation of 
a change in the \t\ dependence of da(pp — > pp)/dt at 
center-of-mass energy y/s — 1.96 TeV. A similar mea- 
surement in pp collisions at y/s — 7.0 TeV has recently 
been reported by the TOTEM Collaboration [13] show- 
ing similar trends for the slope of the differential cross 
section and for the position of the local diffraction min- 
imum, although the local minimum found in pp elastic 
scattering is much more pronounced than the kink we 
observe. 

This article is organized as follows. First we describe, 
in Sec. II, the DO detector, with particular emphasis on 
the FPD. Next, in Sec. Ill, we discuss various aspects of 
the selection of the sample of candidate elastic scatter- 
ing events, including the details on the position measure- 
ments in the FPD, their alignment, the reconstruction 
of p and p scattering angles, the background subtraction, 
the measurements of acceptance and efficiencies. System- 
atic uncertainties are discussed in Sec. IV and the results 
of this measurement are described in Sec. V. 



II. DO DETECTOR 

We briefly describe the elements of the detector that 
are relevant for the measurement reported here. A de- 
tailed description of the DO detector can be found in 
reference [10]. The central tracking system of the DO 
detector comprises a silicon microstrip tracker (SMT) 
and a central fiber tracker (CFT), surrounded by a 2 T 
superconducting solenoidal magnet. The pseudorapid- 
ity [14] coverage for the tracking detectors is \r]\ < 3 
for the SMT and | rj | < 2.5 for the CFT. Outside of the 
superconducting magnet, the liquid- argon and uranium 
calorimeter is composed of three sections housed in sepa- 
rate cryostats: a central calorimeter section covering the 
pseudorapidity range \rj\ < 1.1 and two end calorimeter 
sections that extend coverage to rss 4.2 [15]. A muon 
system [16], outside of the calorimetry, consists of three 
layers of tracking detectors and scintillation trigger coun- 
ters, and large toroidal magnets for muon momentum 
measurement. The luminosity monitor (LM) consists of 
plastic scintillator arrays located at z = ±140 cm (where 
z is measured from the IP along the nominal direction of 
the proton beam), and covers the pseudorapidity range 
2.7 < \rj\ < 4.4. The LM is used to detect non-diffractive 
inelastic collisions and to make an accurate determina- 
tion of the luminosity. 



A. Forward Proton Detectors 

Figure 1 shows the layout of the main components 
of the FPD system relevant to this measurement [10]. 
In the center of the diagram is the interaction point, 
IP, at the center of the main DO detector. A scattered 
proton/antiproton goes through 3 Tevatron quadrupolc 



magnets (with a field gradient of about 20 T/m) which 
alternate defocusing in the horizontal and vertical planes, 
passes through the first station of detectors, goes through 
a region free from magnetic field where electrostatic sep- 
arators are located, and arrives at a second station of 
detectors. 

The FPD consists of four quadrupole spectrometers on 
both the scattered proton (P) and scattered antiproton 
(A) sides plus a dipole spectrometer (not shown). Each 
quadrupole spectrometer is composed of two scintillating 
fiber detectors, one located at about 23 m (Ai or Pi) 
and the second at about 31 m (A2 or P2) from the IP 
along the Tevatron beam line. Both detectors are either 
above (U), below (D), on the inner side (I), or on the 
outer side (O) of the beam line. Given the location of 
the two detectors, the spectrometers on the A side are 
named A\j , Ad , Ai , and Ao (a similar definition is made 
for the P side). The pseudorapidity range covered by the 
detectors is about 7.3 < \r/\ < 8.6. 

Scattered protons and antiprotons cross thin stainless 
steel windows at the entrance and exit of a vessel (Roman 
pot) containing the detectors [17]. The pots are remotely 
controlled and moved to within a few millimeters of the 
beam during stable beam conditions. The Roman pots 
house position detectors using 0.8 mm x 0.8 mm square 
double-clad polystyrene scintillating fibers to detect the 
passage of charged particles. Each position detector con- 
sists of six layers of scintillating fibers (U, U', V, V, X, 
X'), where the scintillating fibers of the U, U' (V, V) lay- 
ers are rotated by plus (minus) 45 degrees with respect 
to the vertical X, X' fibers (see Fig. 2). 

The fibers of the "primed" planes are offset by 0.53 mm 
(two-thirds of a fiber width) with respect to the fibers of 
the "unprimcd" layers. By combining the fiber informa- 
tion from "primed" and "unprimed" layers, we obtain 
"wide" fiber segments (about 1.07 mm wide) that are 
used for triggering and "fine" fiber segments (about 0.27 
mm wide) used for offline hit reconstruction (see Fig. 3). 
Each detector also contains a scintillator (read out by a 
time-to-digital converter system [10]) which provides a 
time measurement for particles passing through the de- 
tector with a resolution of about 1 ns. The time measure- 
ment is used to distinguish particles coming from the cen- 
ter of the DO detector from background beam halo par- 
ticles (particles traveling far enough outside of the main 
beam core that they pass through the FPD detector). 

Beyond the active areas of the detectors, matching 
square clear fibers transport signals from the scintillators 
to 16-channel multi-anode photomultiplicr tubes. The 
electronic signals are subsequently amplified, shaped, and 
sent to the DO triggering and data acquisition system [10]. 

The elastic triggers are defined by the logical OR of 
coincident hits in all four detectors, in either of the four 
possible configurations of a p spectrometer and a diago- 
nally opposite (collinear) p spectrometer: AtjPd, AdPu, 
AjPq, and A Pi- Several different conditions on the 
hits in the scintillating fiber detectors were used in the 
triggers: a tight (T) trigger that registered a single hit 
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FIG. 1: Schematic view of the Roman pot stations (Ai, A2, Pi, P2) comprising the forward proton detector as described in 
the text (not drawn to scale). The dipole spectrometer is not shown. The letters U, D, I, O make reference to the up, down, 
inner and outer detectors, respectively. 




FIG. 2: Schematic view of the design of one FPD scintillating 
fiber detector. U, U', V, V', X, and X' are the scintillating 
fiber layers. The local detector x, y coordinates are indicated 
by the arrows. The scintillating fibers are indicated by the 
green stripes. 

formed by the coincidence of UU', W, XX' wide seg- 
ments; a medium (M) trigger that allowed up to three 
wide segments; and a loose (L) trigger that allowed hits 
formed from coincidences of two out of the three UU', 
VV', XX' wide segments with no requirements on num- 
ber of hits. To reduce backgrounds from inelastic colli- 
sions, LM vetoes (no hits in either the LM counters on 
the proton and antiproton sides) were included as part 
of the elastic triggers. The timing scintillator in each 
detector is only used for providing time information for 
off-line analysis and is not part of the triggers. 



1.07 mm 



0.8 mm 


wide segment 














U fiber layer 



0.53 mm 



0.27 mm fine segment 

FIG. 3: Schematic view of part of the U and U' fiber layers 
with the definitions of wide and fine fiber segments. Similar 
definitions are used for the X, X' and V, V' fiber layers. 

respect to the beam and then use the beam transport 
matrices [18] (which are functions of the currents of the 
magnets located between the IP and the FPD detectors 
and correlate the x, y coordinates and scattering angle 
of a particle at two specific z locations) to reconstruct 
the paths of protons and antiprotons through this re- 
gion of the Tevatron. Next, background subtraction and 
efficiency corrections are performed. We use a Monte 
Carlo (MC) event generator to apply corrections for ac- 
ceptance, detector resolution, beam divergence and IP 
size effects. 



A. Data Sample 



III. ELASTIC EVENT SELECTION 

An initial data sample is obtained by requiring events 
to satisfy one of the elastic triggers. The p and p hit 
coordinates are measured in the FPD system using the 
fibers information and then used to select the sample of 
elastic scattering events. We align the detectors with 



The data for this analysis were collected with dedicated 
beam conditions designed to facilitate the positioning of 
the FPD Roman pots as close to the beam axis as possi- 
ble. The Tevatron injection tune with the betatron func- 
tion of (3* = 1.6 m at the DO IP was used instead of 
the standard (3* = 0.35 m lattice [19]. Additionally, only 
one proton bunch and one antiproton bunch were present 
in the Tevatron. Scraping in the vertical and horizontal 
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planes to remove the halo tails of the bunches was per- 
formed and the electrostatic separators were turned off 
before initiating collisions. The initial instantaneous lu- 
minosity was about 0.5 x 10 30 cm _2 s _1 with a lifetime 
of about 30 hours, corresponding to a mean number of 
w 0.8 interactions per bunch crossing. The recorded lu- 
minosity is about C = 31 nb -1 , which corresponds to the 
sum of two data sets used in this analysis, each with dif- 
ferent detector positions with respect to the beams. One 
set of data was taken with the closest detector position 
reaching about 4 mm with respect to the beam (data 
set 1), and the other data set was taken with detectors 
about 1 mm closer to the beam (data set 2). For the 
given instantaneous luminosity and the conditions of this 
store, about 33% of the elastic events are expected to be 
produced together with an inelastic collision in the same 
bunch crossing. 

Approximately 20 million events were recorded using 
a special trigger list optimized for diffractive physics, in- 
cluding triggers for elastic, single diffractive, and double 
pomeron [20] exchange. This analysis uses elastic trig- 
gers, which make up about 10% of the total data col- 
lected. Independent triggers were used to determine effi- 
ciencies. 



B. Hit Reconstruction 

We reconstruct the proton and antiproton hits using 
fibers with an associated signal above the same thresh- 
old as used in the FPD trigger system, that is tuned to 
accept hits from a minimum ionizing particle. For each 
fiber in an un-primed layer above threshold, we verify if 
there is an adjacent fiber in the primed layer in order to 
define a "fine" fiber segment through which the particle 
traversed the plane. Events with more than four fibers 
firing in a plane are typically not due to single particles 
and are rejected (more than 99% of the clastic triggered 
events survive this condition). We require at least two 
out of the three UU', VV, and XX' fine segments to 
be reconstructed in each detector, with their intersection 
yielding the transverse coordinates of the hit. Since noise 
in some of the fibers could produce fake hits, in the case of 
more than one reconstructed hit in a detector (which hap- 
pens in about 10% of the events), we weight the hits by 
the sum of the analog-to-digital converter (ADC) pulse 
height of all the fibers that contributed to each hit and 
retain only the hit with the highest weight. The cor- 
rection for the selection efficiency, discussed in Sec. Ill I, 
accounts for elastic events discarded by this requirement. 

Utilizing events with fine segments in U, V, and X, we 
are able to estimate the offset in position and resolution 
within each detector by taking the difference between the 
x coordinate obtained from UV fiber intersection and a 
similar measurement using the X fiber plane. Most of the 
detectors used in this analysis have a resolution of about 
150 fxm with the exception of some of the A2 detectors 
that had several inactive scintillating fibers, degrading 
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FIG. 4: Comparison of detector y coordinates in the spec- 
trometers a) Pu and b) Pd- 



their resolution to about 250 fim. 



C. Selection of Candidate Elastic Events 

Elastic collisions produce one of the four possible hit 
configurations in the FPD: AuP D , A D Pu, AjPq, AqPi- 
Since the inner and outer detectors are farther from the 
beam than the U and D detectors, they have poorer ac- 
ceptance for elastic events and are only used for align- 
ment purposes (see Sec. HID). If wc compare the hit 
coordinates reconstructed by each of the two detectors 
of one spectrometer, we observe a correlation band from 
particles going through the spectrometer, but also ob- 
serve some uncorrelated background hits (see Fig. 4). We 
require hits to lie within ±3cr (a rj 220 /zm) from the cen- 
ter of the band. We determine the quantity a by fitting 
a Gaussian distribution obtained by projecting the hits 
onto an axis perpendicular to the correlation band. 

We then investigate the correlation between the co- 
ordinates of the protons and antiprotons in diagonally 
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opposite spectrometers. In addition to the expected cor- 
relation due to the collinearity of elastic events, some 
background contamination due to halo particles remains. 
This contribution is reduced through the use of the FPD 
timing system, but there is still some residual back- 
ground, partially due to inefficiency of the trigger scintil- 
lators and to different acceptances of the spectrometers. 
A correction for the contribution from halo within the 
correlation band is discussed in Sec. Ill F. 

For a specific set of p and p spectrometers in an elastic 
combination, we observe that elastic events with a pro- 
ton passing through a small region of the detectors in one 
spectrometer will have p positions distributed over a simi- 
lar but somewhat larger region of the diagonally opposite 
spectrometer. Depending on the location of the proton, 
the Gaussian distribution of the p position may be trun- 
cated due to the finite detector size. After calculating the 
p acceptance correction as a function of the coordinates 
of the proton spectrometer, a fiducial cut is applied in 
the Pi detectors such that only regions in these detec- 
tors which require a correction of 2% or less are used. 
The region we select in the Pi detectors is thus guaran- 
teed to correspond to an acceptance greater than 98% in 
the other three detectors. Figure 5 shows the y p vs yp co- 
ordinate correlation plots (y coordinate measured in the 
local coordinate system shown in Fig. 6), the tagging 
of the hits according to time-of-flight information, and 
the fiducial cuts applied (indicated by the dashed lines). 
We also check that there is no activity in the calorimeter 
for events in the elastic data sample. The contribution 
to the elastic sample from events that have a non-zero 
energy in the calorimeter is less than 0.1% of the total 
number of selected events. Most of the elastic events pro- 
duced in beam crossings with multiple pp interactions are 
suppressed by the vetoes on the LM detectors. 



D. Forward Proton Detector Alignment 

In order to reconstruct the tracks from the protons 
and antiprotons, we first align the detectors with respect 
to the beam and then determine the position and angle 
of a particle at the IP from the measured hit positions 
in the FPD detectors after application of the Tevatron 
transport matrices. 

The location of the detectors with respect to the beam 
is determined using a sample of tracks that pass through 
one vertical and one horizontal detector at the same Ro- 
man pot station, allowing determination of the relative 
alignment of the detectors. We use elastic events to align 
one horizontal detector that did not have any overlap 
with the other detectors at the same z location. We de- 
fine the center of a vertical and a horizontal pot at each 
Roman pot station and then measure the positions of all 
four detectors with respect to this reference system. Due 
to the beam optics the x (y) hit coordinate distributions 
in the vertical proton (horizontal antiproton) detectors 
are narrower than in the vertical antiproton (horizon- 



tal proton) detectors. Because of this effect, the offset 
of the beam in x for the vertical proton detectors (and 
the offset of the beam in y for the horizontal antiproton 
detectors) is determined by fitting a Gaussian distribu- 
tion to the x and y distributions in these detectors. All 
other beam offsets are obtained given the fact that the 
reconstructed IP offset and scattering angle is the same 
for both the proton and antiproton in an elastic event. 
We take the average over all four elastic combinations 
to determine the offset of the beam at each Roman pot 
station. Figure 6 shows the position of each detector ob- 
tained with the alignment procedure described above for 
the data set corresponding to the detector configuration 
with Roman pots located closest to the beam. The coor- 
dinates are plotted with respect to the beam coordinate 
system (dashed lines shown in the figure). For the data 
collected with the Roman pots retracted further away 
from the beam line, we use positional difference infor- 
mation obtained from the pot motion system added to 
the previously determined aligned positions. The uncer- 
tainty of the location of the detectors after the alignment 
procedure is estimated to be about 200 ^m. 



E. Reconstruction of p/p tracks and |£| 
measurement 

The Tevatron transport matrices are unique for these 
data due to the use of the injection lattice. The beam 
radius (at ler level) ranges from 0.4 mm to 0.8 mm at 
the different FPD detector locations while the beam di- 
vergence is about 44 /xrad. The change in the proton or 
antiproton divergence caused by the quadrupoles system 
is typically of the order of 1 mrad for the smallest \t\ 
values considered in this measurement. 

The size of the IP in x, y and z is typically determined 
by reconstructing the primary vertex using the tracking 
detectors [10] . Given the absence of central tracks in elas- 
tic events, data from non-elastic triggers, simultaneously 
collected with the elastic scattering events used in this 
analysis, are used to obtain the x, y and z distributions 
of the IP. The size of the IP in each direction is obtained 
by fitting a Gaussian function to each distribution. The 
measured Gaussian width of the IP distribution is about 
100 pm in the transverse plane and about 45 cm along 
the z axis (see Fig. 7). The offsets observed in x and 
y from the primary vertex reconstruction are due to the 
fact that the center of the tracker is shifted with respect 
to the beam line. We use the reconstructed elastic events 
to determine these offsets, as described below. 

To tag an elastic event, we reconstruct the coordinates 
of the hits in the two proton and two antiproton detec- 
tors. The difference of coordinates between a proton and 
antiproton detector yields a Gaussian distribution with a 
width related to the fiber offset and resolution of the two 
detectors, the IP size, and the beam divergence. Given 
the previously estimated detector resolutions and IP size, 
these distributions can be used to estimate the beam di- 
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FIG. 5: The correlation plots y p vs pp for the first detectors in the spectrometers (y coordinates measured in the local coordinate 
system shown in Fig. 6) for AuPd and AdPu- Plots a) and d) show the correlations without any timing requirement; b) and 
e) correspond to hits with a coincident tag in the corresponding scintillator early time windows; c) and f) show events with no 
hits within the early time window (most of the hits in this case are in the in-time window) . The dashed lines correspond to 
the fiducial requirements applied. 
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FIG. 7: Size of the IP determined by the reconstruction of the 
primary vertex from non-elastic triggers (black points) fitted 
with a Gaussian function (solid line), a) x distribution of the 
IP. b) y distribution of the IP. c) z distribution of the IP. 



FIG. 6: Detector positions with respect to beam center 
(dashed lines) for the data set corresponding to the closest 
pot insertion to the beam. The arrows indicate the local co- 
ordinate system in each detector. 

vergence. The values obtained are similar (40 ± 5 /irad) 
to those estimated using the injection tune lattice param- 
eters. With the coordinates of the two detectors in each 
spectrometer we determine the offsets of the IP (xq, yo) 



and the horizontal and vertical scattering angles (9 X , 9 y 
respectively) by using the following equations: 

Xi = M Xji x + L Xti 9 x (1) 

where Xi, yi are the hit coordinates in a detector i and 
M Xt i, L x i (My i, L y i) are the transport matrix elements 
in the horizontal (vertical) axis for that detector loca- 
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FIG. 8: Difference in the reconstructed \t\ between p and p 
vertical spectrometers. The MC (solid line) expected distri- 
bution is compared to data (solid points). a)AuPo combina- 
tion, b) AdPu combination. The number of entries in each 
distribution is normalized to unity. 

tion. We verify that the distributions of the difference 
of the IP offsets and scattering angles obtained from the 
proton spectrometer and the antiproton spectrometer for 
every event are centered around zero. Since the values 
of 8 obtained are of the order of milliradians, the four- 
momentum-transfer squared, t, can be approximated as 

t=- P \el + el) (2) 

where p is the momentum of the scattered particle. Given 
that the momenta of the elastically scattered proton and 
antiproton are the same as that of the incoming parti- 
cles (980 GeV) and that the beam momentum spread is 
negligibly small (0.014%) [21], the uncertainty in \t\ is 
dominated by the measurement uncertainty of the scat- 
tering angle. Figure 8 shows the difference in the re- 
constructed \t\ from the p and p vertical spectrometers, 
A|i| = \t\p — \t\ p . The non-zero standard deviation of 
cr&t, is due to the resolution on \t\ p and \t\p. The 
average of \t\ p and \t\ p (\t\ avc = (\t\ p + \t\ p )/2) has a reso- 
lution of approximately cta|*|/2. The \t\ bin size is chosen 
as the largest of the observed &A\t\ for the two detector 
combinations AtjPd and AdPu, where the difference in 
<?A\t\ results from different detector resolutions. We also 
study the \t\ resolution as a function of \t\ and observe a 
gradual increase with from 0.02 GeV 2 to 0.04 GeV 2 , 
which is reflected in the bin widths. 



F. Background Subtraction 

The primary source of background in the selected sam- 
ple is due to beam halo, consisting of either a halo pro- 
ton and a halo antiproton in the same bunch crossing or 
a halo particle combined with a single diffractive event. 
A halo particle passing through the proton detectors in 



the time window for the protons which have undergone 
elastic scattering usually passes through the diagonally 
opposite antiproton detector at an earlier time (and vice 
versa). Therefore, the time information can be used to 
veto events with early time hits, consistent with halo pro- 
tons and antiprotons in the elastic sample. The veto is 
not 100% efficient due to a combination of scintillator 
efficiency and positioning of the detectors with respect 
to the beam (a closer detector position both increases 
a detector's signal acceptance and its ability to reject 
halo). Consequently, it is necessary to subtract the re- 
maining background. We consider an event to be caused 
by P (p) halo if one or both of the two p (p) detectors 
of an elastic combination have hits in their early time 
interval. We require events to have no activity in either 
the p ha i or the Phaio timing window. We select back- 
ground samples by requiring hits consistent with p halo 
and p halo simultaneously. First, we verify that outside 
the elastic correlation band between the coordinates of a 
proton and antiproton detector the signal tagged events 
have the same \t\ dependence as the background tagged 
ones (see Fig. 5). Next, assuming that signal and back- 
ground tagged events also have the same \t\ dependence 
inside the correlation band, we use the ratio of these two 
distributions and use it to estimate the percentage of 
background events inside the signal tagged correlation 
band. This background is subtracted from all events in- 
side the correlation band to obtain dN jdt as a function 
of |t|. The amount of background subtracted inside the 
elastic correlation band varies from 1% at low \t\ to 5% 
at high \t\. The absolute uncertainty of the background, 
which is propagated as a statistical uncertainty to da/dt, 
varies from 0.3% at low \t\ to 5.0% at high \t\. As a cross 
check, we vary the detector band cuts from 3. Oct to 3.5ct 
and to 6. Oct, to allow more background, and obtain sim- 
ilar dN/dt results after applying the same background 
subtraction procedure (within 1%). 



G. Monte Carlo Simulation of Elastic Events 

We have developed a MC generator interfaced with the 
Tevatron transport matrices to generate elastic events. 
The MC allows us to study the geometrical acceptance 
of the detectors, resolution of the position measurement, 
alignment, and effects of the beam size and beam diver- 
gence at the IP. The generation of events is based on an 
Ansatz function that we obtain by fitting the dN/dt dis- 
tribution of the data. We study the acceptance and bin 
migration effects using samples generated with a wide 
range of different dN/dt Ansatz distributions. The vari- 
ations in the corrections are included as systematic un- 
certainties. 

The positions of hits in each detector are transformed 
into fiber hit information, and the reconstruction then 
proceeds using these hits, following the same procedure 
as with the data. The reconstructed correlation patterns 
in MC are in good agreement with those observed in data. 
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FIG. 9: Comparison of MC (solid line) and data (black points) 
for the coordinate difference y P - for the elastic configura- 
tion Ai,tjPi,d for the data set with detectors closer to the 
beam. 
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FIG. 10: Azimuthal acceptance of the FPD in data, after 
fiducial cuts, for the closest detector position of the AuPd 
elastic combination. The points correspond to MC, the solid 
line corresponds to the calculation of the <f> acceptance from 
the detector positions and the fiducial area. 



In addition, the MC also predicts the widths of the dif- 
ferent correlations as shown in Fig. 8 and Fig. 9 . 



I. Selection and Trigger Efficiencies 



H. Acceptance and Bin Migration Correction 

Due to accelerator optics, protons and antiprotons of a 
particular \t\ are transported to an elliptical region in the 
x — y plane at each detector location. Each detector has a 
different coverage in the azimuthal angle (j> and therefore 
a different coverage of the ellipses of constant \t\. The 
4> acceptance correction accounts for the fraction of the 
\t\ ellipse for each \t\ bin that is not covered by the fidu- 
cial area used in each detector. This correction depends 
only on the detector location with respect to the particle 
beam and on the geometry of the fiducial area used in the 
detector for selecting the events. Figure 10 shows the <j) 
acceptance calculated as a function of |i| for data at the 
closest detector position with respect to the beam and 
after the fiducial cuts. We obtain similar results for the 
<f> acceptance using the MC described in the previous sec- 
tion before adding the effects of beam divergence, IP size 
and detector resolution. The uncertainty in the <p accep- 
tance correction, which comes from the size of the MC 
sample used, is less than 0.1% at low \t\ values and less 
than 1.0% at high \t\ values. To estimate the correction 
for bin migration effects that is applied in data to obtain 
the da/dt distribution, we include the measured values 
of beam divergence, IP size, detector resolution and pot 
position uncertainties in the MC generator and compare 
to a same size MC sample with no smearing effects. 



We simultaneously determine the effects of selection 
and trigger efficiencies for each of the four detectors cor- 
responding to an elastic configuration. To determine the 
efficiency of a particular detector, we use an indepen- 
dent trigger which does not include that detector. We 
obtain the dN/dt distributions with an elastic track re- 
constructed in the other three detectors and for all four 
detectors, with the trigger conditions satisfied. The ratio 
of the two distributions is used to extract the efficiency of 
the detector as a function of \t\ (where \t\ is reconstructed 
from the coordinates of the opposite side spectrometer). 
To avoid any additional effect from detector acceptance, 
we use only hits in the other three detectors in a region 
determined to have full geometrical acceptance in the 
detector of interest. We repeat a similar procedure for 
each of the four detectors in every elastic combination 
and multiply the efficiencies of the four detectors to de- 
termine the final efficiency correction. Typical selection 
and trigger efficiencies are in the range of 50% to 70% 
depending on the detector and trigger requirements. Ad- 
ditionally, we make a correction for the veto in the LM 
which was part of the elastic trigger. This veto filters 
out elastic events that are produced in coincidence with 
an inelastic collision in the same bunch crossing (pileup). 
To make this correction, we use a trigger based on hits in 
the proton spectrometers with no LM requirement and 
determine the fraction of candidate elastic events recon- 
structed in coincidence with the LM. We find that about 
27% (20%) of elastic events for data set 1 (2) as defined 
in Sec. Ill A were removed by the LM veto. 
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J. Luminosity 

Since elastic data were collected with Tevatron con- 
ditions modified with respect to standard operations of 
the DO experiment, we compare the number of inclusive 
jet events obtained from our data to the number of jet 
events used for the inclusive cross section measurement 
discussed in [22] to determine the integrated luminosity. 
The corresponding integrated luminosity for data set 1 
(2) is 18.3 nb" 1 (12.6 nb" 1 ), with an uncertainty of 13%. 
We add in quadrature the uncertainty in the standard lu- 
minosity determination (6.1% [23]) and obtain an overall 
normalization uncertainty of 14.4%. 

IV. SYSTEMATIC UNCERTAINTIES 

The major |t|-dependent contributions to the system- 
atic uncertainties in the measurement of da/dt are due 
to detector efficiencies, beam divergence, detector posi- 
tions, and the choice of the Ansatz function. The lumi- 
nosity measurement contributes to the overall normaliza- 
tion uncertainty. The Tevatron beam transport matrices 
are known with high precision (within 0.1%) and there- 
fore produce a small uncertainty in our results compared 
to the other sources. We take the uncertainty in the posi- 
tion of the pot from the beam center as an extra "smear- 
ing" factor for the hit coordinates in the MC. Since the 
efficiencies vary with |f|, we fit either a polynomial or an 
exponential function to each trigger efficiency and prop- 
agate the uncertainties in the fit parameters to da/dt, 
using the covariance matrix of the fit. For the beam di- 
vergence term, we vary the beam divergence by ±5 /irad 
in the MC, and we propagate the change in the accep- 
tance correction bin-by-bin to da/dt. We also consider 
26 possible variations of the Ansatz function used in the 
MC to account for the uncertainties in the logarithmic 
slopes before and after the kink and also for the \t\ value 
where the kink is observed. 



V. RESULTS 

In total, we have four independent measurements 
of da/dt that come from the two elastic combinations 
(AuPd and AdPu) and two data sets, which agree with 
each other within the uncertainties. We combine the four 
measurements using a bin-by-bin weighted average. The 
resulting values for the da/dt distribution, together with 
their total uncertainties, are listed in Table I and shown 
in Fig. 11. The uncertainties are the total experimen- 
tal uncertainties excluding the 14.4% normalization un- 
certainty The \t\ bin centers are determined using the 
prescription described in [24], however, the values found 
are very close to the middle of the bin. Two phenomeno- 
logical model predictions for y/s = 1.96 TeV (BSW [25], 
Islam et al. [26]) are also shown. The BSW model shows a 
good description of the data in shape and normalization 



TABLE I: The da/dt differential cross section. The statistical 
and systematic uncertainties are added in quadrature. The 
luminosity uncertainty of 14.4% is not included. 
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FIG. 11: The measured da/dt differential cross section. The 
normalization uncertainty of 14.4% is not shown. The uncer- 
tainties are obtained by adding in quadrature statistical and 
systematic uncertainties. The predictions of BSW ([25]) and 
Islam et al. ([26]) are compared to the data. 



and is able to reproduce the kink within experimental 
uncertainties. 

The \t\ range covered by our measurement is 
0.26 < |t| < 1.2 GeV 2 . We observe a change in 
the logarithmic slope of the da/dt distribution at 
\t\ « 0.6 GeV 2 . A fit to the da/dt distribution in 
the range 0.26 < \t\ < 0.6 with an exponential func- 
tion of the form Ae~ b ^ yields a logarithmic slope pa- 
rameter of b = 16.86 ± 0.10 (stat) ± 0.20 (syst) GeV" 2 
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FIG. 12: The do/dt differential cross section measured by the 
DO Collaboration and compared to the CDF and E710 mea- 
surements at y/s = 1.8 TeV, and to the UA4 measurement at 
y/s = 0.546 TeV (scaled by a factor of 100). A normalization 
uncertainty of 14.4% on the DO measurement is not shown. 

(x 2 = 6.63 for 7 degrees of freedom). 

Figure 12 shows a comparison of our results to those 
obtained at y/s = 1.8 TeV by the CDF and E710 Teva- 
tron Collaborations [11, 12], and also to that of the UA4 
Collaboration at y/s = 0.546 TeV [27]. Our measurement 
of the slope parameter agrees within uncertainties with 
previous measurements by the CDF (b = 16.98 ± 0.25 
GeV" 2 ) and E710 (b = 16.30 ± 0.30 GeV" 2 ) Collabora- 



tions. A comparison of the shape of our measured da/dt 
to UA4 measurement shows that the kink in da/dt moves 
towards lower |i| values as the energy is increased, but, 
as in the UA4 data, we do not see a distinct minimum as 
observed in pp elastic interactions ([13]). 

In summary, we have presented the first measurement 
of da(pp — > pp)/dt as a function of \t\ at y/s = 1.96 TeV. 
Our measurement extends the |t| and y/s range previ- 
ously studied and shows a change in the |t| dependence, 
consistent with the features expected for the transition 
between two diffractive regimes. 
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